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We show that a physical realization of the phase diagram proposed by Minnhagen et al. [Phys. Rev. B 78, 
184432 (2008)] for so-called generalized fully frustrated XY model on square lattice can be achieved in arrays 
of SFS (superconductor-ferromagnet-superconductor) junctions near the transition of the junctions to ;r-state. 
Moreover, the phase diagram with such a topology has to be reproduced twice, on both sides of the 0 - ;r 
transition. 


I. INTRODUCTION 

A fully frustrated (FF) XY model can be defined by the 
Hamiltonian 

(la) 

UJ') 

where the form of the interaction of variables 6j dehned on 
sites j of some regular two-dimensional lattice is described 
by an even periodic function U{(f>) (with period In) minimal 
at (^ = 0. The summation in Eq. (la) is performed over all pairs 
(jf) of nearest neighbors on the lattice. The non-fluctuating 
(quenched) variables s -Ayj defined on lattice bonds 
have to satisfy the constraint 

- +7T(modln) (lb) 

□ 

on all lattice plaquettes. The notation □ below the sign of 
summation implies the directed sum of variables Ajy = -Aj>j 
over the perimeter of a plaquette in the counterclockwise di¬ 
rection. The model is gauge-invariant, that is, its properties do 
not depend on the particular choice of gauge variables Ajf as 
long as constraint (lb) is satisfied. 

The model defined by Eqs. (1) with interaction function 

t/(«^) = t/o(<^) = 7(l-cos<^) (2) 

can be used for the description of a planar magnet with odd 
number of antiferromagnetic bonds per plaquette [ 1 ] and of a 
magnetically frustrated Josephson junction array with a half¬ 
integer number of flux quanta per plaquette [2]. In the lat¬ 
ter case, 6j is the phase of the order parameter on the y'th 
superconducting island and (pjf = Of - Oj - Ajf the gauge- 
invariant phase difference between the neighboring islands. 
Note that here and below Uitp) is always assumed to be 
counted off from its value at 0 = 0, where U((p) is minimal. 

In contrast to the case of the FF XY model on dice lattice 
[3], the ground states of the FF XY model on square or tri¬ 
angular lattice can be constructed by minimizing the energy 
of each plaquette independently. For a square plaquette, the 
minimization of 

4 

Fplaq - (3) 

a=I 

under the constraint - am following from Eq. (lb) 

for t/(0) gives 

(4) 


which immediately allows one to construct a ground state with 
a checkerboard structure which in addition to continuous U(l) 
degeneracy (related to the simultaneous rotation of all vari¬ 
ables 0j) has also a two-fold discrete degeneracy related to the 
choice of sign in Eq. (4) [1]. 

In accordance with that, the model allows for the existence 
of two phase transitions related to the disordering of the con¬ 
tinuous and discrete degrees of freedom (see Refs. 4 and 5 for 
reviews). The former takes place at slightly lower temperature 
than the latter due to the screening of the interaction of the 
conventional (integer) vortices by the kinks on domain walls 
behaving themselves as fractional vortices [6]. 

Few years ago Minnhagen et al. [7, 8] suggested that a sub¬ 
stantial change in the form of the interaction function in Eq. 
(la) may lead to a change in the sequence and/or nature of 
phase transitions in the system. They chose to replace Uoitp) 
by the interaction function (5) introduced earlier [9, 10] for 
studying the transformation of the phase transition in the con¬ 
ventional XY model from a continuous to a first-order one. 
The phase diagram of the FF XY model with such an inter¬ 
action (named the generalized FF XY model) has been in¬ 
vestigated in Refs. 7 and 8 by numerical methods. It has a 
rich structure with four different phases separated by transi¬ 
tion lines of various types and incorporates four multicritical 
points. 

The main aim of the present note is to point out that a phys¬ 
ical realization of the phase diagram proposed by Minnhagen 
et al. for the generalized FF XY model on square lattice can 
be achieved in arrays of SFS (superconductor-ferromagnet- 
superconductor) junctions [11-13] near the transition of the 
junctions to so-called 7r-state [14-16]. After analyzing in Sec. 
II the main reasons for the overall structure of the phase dia¬ 
gram of the generalized FF XY model, in Sec. Ill we demon¬ 
strate that in the most relevant range of parameters the func¬ 
tion describing the energy of a SFS junction has almost the 
same form as the one used in Refs. 7 and 8. As a conse¬ 
quence, a magnetically frustrated array of SFS junctions with 
half-integer number of flux quanta per plaquette has to demon¬ 
strate a rich phase diagram with the structure analogous to that 
of the generalized FF XY model. 


II. GENERALIZED EE XY MODEL 

As has been already mentioned in the Introduction, in Refs. 
7 and 8 Minnhagen et al. investigated the FF XY model de- 
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fined by Eqs. (1) with interaction function. 
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This interaction function had been introduced in Ref. [9] to 
study the influence of the form of the interaction (depending 
on the value of parameter p) on the type of the phase transition 
in the conventional (without frustration) XY model. At p = 1 
function t/gen(</’) defined by Eq. (5) coinsides with Uo((f>). The 
p-dependent numerical prefactor in the right-hand side of Eq. 
(5) had been introduced to make t/"(0) = J independent of p. 
This parameter determines the energies of smooth continuous 
fluctuations (spin waves and vortices) in the unfrustrated case. 
The EE XY model with interaction function (5) was named in 
Ref. 7 the generalized fully frustrated XY model. 

The most important property of the generalized EE XY 
model on square lattice (responsible for the overall structure 
of its phase diagram) is the change in the structure of its 
ground states [7], which takes place with the increase in p 
at 




Pc = 


ln(3/4) 


I 21n[cos(;r/8)] 


1.347880 . 


( 6 ) 


At p < Pc the ground states of the generalized EE XY model 
have the same structure as in the conventional one (corre¬ 
sponding to p = 1). However, at p > pc the energy of a single 




FIG. 1. (color online) (a) Comparison of the form of t/gen(0) for p=l 
(upper curve), p = 1.2 and p = iA (lower curve); (b) comparison of 
the form of t/sFs(<^) for t = 0 (upper curve), t = 0.125 and t = 0.25 
(lower curve). Vertical dashed lines correspond to 0 = ;r/4 and <p = n. 
In this and all other plots the energies of a bond (or of a plaquette) 
are given in units of coupling constant J. 


FIG. 2. The dependence of plaquette energy i?piaq(0) on the value of 
gauge-invariant phase difference (p on one of its bonds (a) for the case 
of the generalzied interaction Ugcn(<f>) with p = Pc', (b) for the case of 
two-harmonics interaction Usfs(P) with t = tc- 


plaquette is minimized when on one bond (p - n, whereas on 
three others bonds p -0. 

In order to explain the reasons for the modification of the 
ground state structure in the generalized EE XY model with 
the increase of p, Eig. 1(a) compares the form of t/gen(<^) for 
three different values of parameter p {p - 1, p - 1.2 and 
p = 1.4). It is evident from this plot that the increase of p 
suppresses the value of t/gen(^) much stronger than that of 
t/gen(7r/4). This suggests that instead of uniformly distribut¬ 
ing the gauge-invariant phase difference ±n between the four 
bonds belonging to a plaquette it can become more advan¬ 
tageous to concentrate it completely on one bond with three 
others having the minimal energy, that is 0 = 0. 

Under the condition that the value of gauge-invariant phase 
difference p on one of the bonds differs from that on three 
other bonds (for which the values of p are assumed to be 
equal) the energy of a single square plaquette is given by 

£plaq(<^) = t/(<^) + 3t/[(7r-^i)/3]. (7) 

Eor p < Pc function £’piaq(0) is minimal at 0 = njA (when 
(p - 7r/4 on all four bonds), whereas for p > pc it is minimal 
at 0 = TT (when on three other bonds (p - 0). Ai p - pc, 
the two minima have equal depths, see Eig. 2(a). At this point 
the first-order phase transition between the two families of the 
ground states takes place. Accordingly, the twofold discrete 
degeneracy of the ground states (existing at p < p^), at p > p^ 
is replaced by a much more developed discrete degeneracy 
which leads to a change in the classification of topological 
excitations (and in the nature of low temperature phase) and 
to the emergence of a different scenario for the disordering of 
the system. 
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III. FF XY MODEL WITH TWO-HARMONICS 
INTERACTION 

Let us now check how the form of the interaction function 
evolves if instead of smoothly changing the exponent in Eq. 
(5) one adds to t/o(0) the second harmonics, 

t/(0) = 7i(l - cos (^)-H 72(1 - cos20). (8) 

The two-harmonics function (8) can be used for the descrip¬ 
tion of SFS junctions in the vicinity of so-called 0 - n tran¬ 
sition, at which the amplitude of the first harmonics changes 
sign [14-16]. In this region of parameters, the amplitude of 
the main contribution to the interaction function [having the 
form t/o(0)] is strongly suppressed and may be of the same 
order as that of the next (second) harmonics. Investigation of 
the conventional (without frustration) XY models with two- 
harmonics interaction of the form (8) started long ago [17, 18] 
and became rather active in recent years, see Refs. 19-21 and 
references therein. 

For further analysis it is convenient to rewrite Eq. (8) with 
7 i_ 2 > 0 as 

UsFs(<f>) = ~ ^) + - cos 2(p)] , (9) 

where f = J 2 IJ 1 denotes the ratio of the two amplitudes, and 
J - J\ +AJ 2 is the effective coupling constant which is used in 
the figures as the unit of energy. Like in the case of t/gen(0), we 
have included into the definition of t/sFs(0) a dimensionless 
prefactor whose value is chosen to make - J inde¬ 

pendent of parameter t responsible for the shape of UsfsW- 

Fig. 1(b) demonstrates how the increase of the ratio of the 
amplitudes of the second and first harmonics t changes the 
form of function Usps((p)- Comparison of Fig. 1(b) with Fig. 
1(a) shows that the increase of parameter t influences the form 
of the interaction function basically in the same way as the in¬ 
crease of parameter p in Eq. (5). In particular, the value of 
UsFsi^) is suppressed with the increase in t much stronger 
than that of t/sFs(^/4). This suggests that the FF XY model 
with interaction (9) can be expected to demonstrate the recon¬ 
struction of the ground state structure completely analogous 
to the one in the generalized FF XY model of Ref 7. Indeed, 
this reconstruction takes place at 

tc 0.207108, 



FIG. 3. (color online) Comparison of the form of (7gcn(0) at p = 
(upper curve) and that of Usfs(4’) for ^ - h (lower curve). 



FIG. 4. The form of Ugtn{<f>) at p = 30. Note an extended scale on the 
horizontal axis. 


where the two minima of the function Fpiaq)*^) defined by Eq. 
(7) have equal depths, see Fig. 2(b). Like in the case of t/gen(0) 
with p > Pc, at t > tc the energy of a plaquette is minimal 
when on one of its bonds (ft - n and on the other three bonds 

= 0 . 

To further emphasize the similarity of the two models, in 
Fig. 3 we compare the form of Ugcni<f>) P - Pc that 
of UsFs{<P) at f = tc- This graph demonstrates that these two 
functions are very close to each other. This is not a mere coin¬ 
cidence. Although for an arbitrary p function UgenW defined 
by Eq. (5) looks like an unphysical one, at 

p = p. = V2;^ 1.414214 (10) 

it exactly coincides with function t/sFs(0) defined by Eq. (9) 
with 

f = L=0.25. (11) 

Since pc is relatively close to p*, and L to f*, 

p» - Pc ~ 0.160(p. - 1), L - L ~ 0.172f«, (12) 

it is no wonder that the two functions plotted in Fig. 3 are very 
close to each other. 

At this point it has to be emphasized that the most important 
features of the phase diagram of the generalized FF XY model 
(namely, the line separating the two low-temperature phases 
from each other and three multicritical points) are restricted to 
the narrow (< 1%) interval around pc, see Fig. 5 of Ref. 8. The 
extreme closeness of the two interaction functions depicted in 
Fig. 3 allows one to expect the existence of analogous features 
in the phase diagram of the magnetically frustrated SFS array 
in a narrow vicinity of the line 

4.828. (13) 

72(7’) tc 

Since the two models have the same classification of topo¬ 
logical excitations in both of their low-temperature phases 
(below and above pc or L), the general topology of their phase 
diagrams can also be expected to be the same. However, one 
of the features of the phase diagram of the generalized FF XY 
model, namely, the transformation of the phase transition from 
a continuous to a first-order one taking place at p a; 30 [8], 
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FIG. 5. (color online) The schematic structure of the phase diagram 
of a fully frustrated array of SFS junctions with two-harmonics in¬ 
teraction (8). The phase denoted Z2XU(1) combines a spontaneously 
broken Z 2 symmetry with a finite value of helicity modulus, whereas 
in the phase denoted Z 2 the helicity modulus is equal to zero. The 
phase with a modified structure of the ground states (leading to their 
extensive degeneracy) is denoted N. 


cannot be reproduced in the model with two-harmonics inter¬ 
action (9) because the form of Ugen{<t>) at p 30 is extremely 
flat outside of a very narrow vicinity of (/i = 0 (see Fig. 4) and 
in no way can be approximated by the two-harmonics function 
t^SFs(0)- 


IV. CONCLUSION 

To conclude, in the present note we have shown that a mag¬ 
netically frustrated square array of SFS junctions [11-13] with 
half-integer number of flux quanta per plaquette has to demon¬ 
strate a rich phase diagram with the structure quite analogous 
to that of the generalized FF XY model [7,8] with interaction 
(5). Like in the case of the unfrustrated SFS array with square 


lattice [22], the phase diagram of the fully frustrated square 
SFS array has to be symmetric with respect to the change of 
the sign of parameter Ji in Eq. (8) at constant value of 72- 
This follows from the possibility to shift in Hamiltonian (la) 
all variables Ajy by ji without violating the constraint (lb). 
Accordingly, in the phase diagram of the FF square SFS array 
the phase diagram of the generalized FF XY model has to be 
reproduced twice, on both sides of the 0- ji transition. 

The schematic structure of this phase diagram is shown in 
Fig. 5, where it is plotted in variables J\IJ 2 and 7’/72- In ad¬ 
dition to the three phases which exist in the conventional fully 
frustrated XY model [5] (corresponding to ^2 —> 0 limit) this 
phase diagram includes also an additional phase denoted N. 
Although in Ref. 8 the properties of this phase are character¬ 
ized with the help of a rather complicated order parameter (in¬ 
volving counting the number of comers on domain walls) we 
expect that it can be more transparently described by compar¬ 
ing the behavior of correlation functions for variables exp W 
and exp 2i0, which is relegated to a separate publication. 

A fmstrated XY model with square lattice and two- 
harmonics interaction has been studied earlier in Ref. 23. 
However, for some unexplained reason, the authors of this 
work have chosen to consider an interaction which has not 
gauge-invariant form, which raises doubts on applicability of 
their results for describing physical systems. The model suit¬ 
able for the description of FF SFS array with triangular lattice 
has been studied in Ref. 24, however this work considers the 
case which in the notation of our Eq. (8) corresponds to ^2 < 0 
and cannot be compared to the case J 2 > 0 considered here. 
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